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1. INTRODUCTION 

Distributed generation like solar, wind, mini hydro, and hybrids are gaining popularity in remote 
locations owing to the fact that it needs no transmission lines, is less polluting, and is efficient and economic. 
A solar photovoltaic (SPV) system provides clean power to consumers at a low rate and moderate efficiency. 
With the increase in connection of SPV systems into the grid, it also gives rise to certain limitations, like 
steady state and temporary overvoltage, power fluctuations, harmonics, and other power quality (PQ) issues. 
These PQ issues have been taken care of by power electronic converters such as distribution static 
compensator (DSTATCOM). The DSTATCOM acts as a shunt connected active power filter, which provides 
reactive power compensation as well as harmonic reduction [1]. The PV inverter itself will act as a 
DSTATCOM and hence increase the utilisation of the switches by two fold, which otherwise would remain 
idle in the absence of PV power. Several control strategies for reference current generation and configuration 
of DSTATCOM are reported in various literature. Various maximum power point tracker (MPPT) algorithms 
for maximum power extraction that can be possible from PV arrays under variable solar insolation are also 
reported in [2]-[4]. Adaptive algorithms based on adaptive linear neuron (ADALINE) least means square 
(LMS) and variations on LMS, such as least means forth (LMF), combined LMS-LMF, recursive LMS, and 
VLG-LMS, have gained prominence over conventional control algorithms such as synchronous reference 
frame theory (SRF), instantaneous reactive power theory (IRPT), and SCT as they provide fast response 
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under variable input conditions such as solar and wind. Godavarti and Hero [5] presents the technique for 
updation of the parameters in the LMS filter in a manner that reduces the computational burden. 
Kohli and Mehra [6] discussed the analysis of MG-LMS two-step algorithms under time-varying conditions. 
The setting of smoothing and control parameters improves tracking performance and stability. Authors have 
shown a reduced complexity filtered-s LMS (FSLMS) algorithm for controlling non-linear noise [7]. 
References in [8]-[10] used a different variant of LMS algorithm in which the regularisation parameter is 
updated dynamically and in [11], fast fourier transform (FFT) and fast hartley transform (FHT) is employed 
for reducing computational complexity. Giri et al. [12]-[14] demonstrate variation in the basic LMS 
algorithm to increase system dynamic response and to adapt to variable input conditions. Research by 
Agarwal et al. [15], the learning loss mitigation funding (LLMF) theory is proposed in which a leakage factor 
is introduced to alleviate stalling and stability system response. Khan et al. [16] proposed the least-mean- 
sixth algorithm for PQ mitigation. The Kalman filter approach is discussed in [17]. Harris hawk optimization 
technique is proposed in [18]. The four-wire multilevel configuration is discussed in [19]. Performance of 
QLMS approach for abnormal condition is discussed in [20]. In the presented work, a zero attractor is 
included for updating the tap weight vector, and a lo norm driven tap vector constraint is incorporated in cost 
function of LMS, thus enhancing the convergence speed. The adaptive zero attractor is selected in lo norm 
constraint LMS (lo-LMS) algorithm and the optimal value for each iteration is derived by reducing transient 
mean square deviation. The performance of the PV-DSTATCOM system is observed for active power 
transfer from PV to grid, reactive power as required by load, mitigating PQ challenges like harmonics, unity 
power factor operation, and voltage control for nonlinear and unbalanced load conditions. This algorithm is 
established in the signal-processing domain, and the results of its implementation for addressing PQ issues 
are found satisfactory. The adaptive lo-LMS algorithm is utilised here to solve PQ concerns such as low 
power factor and harmonics reduction. The peak power from the PV module is extracted using the 
incremental conductance (InC) MPPT technology, and the active power is given to the load connected to the 
point of common interface (PCI) as well as into the grid. The adaptive lo-LMS algorithm is used to operate 
voltage source inverter (VSI) for reactive power compensation and to eradicate harmonics from the system. 
Also, the system stability is maintained during unbalanced loading conditions. In this paper, section 2 
contains the system configuration, a description of control scheme, and different modes of operation. The 
simulation results of proposed system are reported in section 3. Section 4 discusses the conclusion. 


2. METHOD 

Figure | depicts the proposed system configuration. It comprises a PV array, a three-phase VSI, an 
interfacing inductor, an RC filter, and a nonlinear load connected at PCC. A single-stage converter topology 
has been chosen to connect the PV-DSTATCOM system to three-phase grid at PCC. Several MPPT 
techniques have already been reported in several studies for obtaiing peak power from a solar PV array under 
diverse environmental circumstances [21]-[23]. In comparison to other techniques, an InC algorithm is found 
to be faster. The peak power from the PV array is extracted by implementing an InC based approach over 
varying environmental circumstances. 
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2.1. Proposed lo-LMS algorithm 

In this work lo norm constraint LMS algorithm with adaptive zero attractor is presented. Instantaneous 
values are taken into account in the standard LMS algorithm that might cause slower convergence. In the 
proposed algorithm, a zero attractor (ZA) term is added, which accelerates the convergence speed by attracting 
near zero coefficient to zero. However, the performance of lo norm is dependent on value of zero attractor term. 
The setting of a fixed zero attractor is difficult when the system has a time-varying noise signal. We present an 
adaptive algorithm to select ZA in the lo norm. Estimation of weights of real and reactive load current 
components, real and reactive loss components are the main functions of the AZA lo-LMS control algorithm. In 
addition, templates for in-phase and quadrature units are extracted and reference grid currents are generated. 
Actual current is then compared with reference in a hysteresis current controller and corresponding gating 
signals are generated for switching the VSI. Consider the following linear regression model. 


d(k) = w' x(k) + v(k) (1) 


Where d(k) symbolized desired signal, w is unknown impulse response vector, x(k) indicates input vector 
whereas noise signal is denoted as v(k). In the proposed lo-LMS algorithm, the weights are updated as, 


wk + 1) = w(k) + pe(k)x(k) + po(w(k)) (2) 


Where w(k)=[Wo(k), wi(k), ..., Wi-1(k)]? denotes filter tap weight vector, k represents the time scale, the parameter 
u is LMS convergence coefficient or step size of adaptation. The error signal e(k) is calculated as in (3), 


e(k) = d(k) — w"(k)x(k) (3) 


The term p&o(w(k)) in (2) represents the zero-attraction term, in which the parameter p > 0 denotes 
the zero attractor, whereas the term €o(w(k)) implies the subgradient of the tap weight penalty in lo norm, and 
the subgradient is specified as in (4) for each entry of 0<=i<L, 


&)(w,(k)) = a?w;(k) — asgn(w;(k)), < < (a) 
= 0, otherwise 


Where the strength as well as range of zero-attraction is governed by positive parameter o>0. With reference 
to steady state mean square deviation (MSD), the optimum zero attractor is positively connected to noise 
power (17) in [24]). If there is any change in the level of the signal, the zero attractor resets to improve the 
steady-state performance of the lo-LMS algorithm. We introduce the adaptive zero attractor along with 
lo-LMS algorithm (AZA lo-LMS) to improve the ability of the algorithm to adapt to time-varying signals. 
Based on the objective of minimising transient MSD as much as achievable during the convergence process, 
we begin by optimising the drop in transient MSD to get the optimum value of zero attractor after every 
iteration. Then practical values are obtained from the optimal solution. The optimal value of ZA is given as, 


— Elw' e)go(w(k))|-mEle(e)x? (K)Fo(w))] (5) 
Pape E[E} w(k)Eo(w(k))] 


The time averages of all three expectation on RHS of (5) is approximated as, 


Wk + 1) = BY) + (1 Bow" (kg (wk) (6) 
b(k + 1) = BOK) + (1 - B)x"(k)So(w(k)) (7) 
P(k +1) = BAK) + (1 BEF w(k)Eo(w(k)) (8) 


Where w(k+1), o(k+1) and o(k+1) represents the three expectation terms on RHS of (5). Furthermore, the 
value of smoothing parameter B is chosen to set to near unity for proper time smoothing. In this paper, the 
value of B is selected at 0.99. The reference optimal ZA is now written as, 


p(kt+1)-ub(k+1) 
pkey k+1) #0 

Popt (k +1)= (k+1) if &( ) - 
0 


if o(k+1) =0 


In case of denominator equals to zero in (9) when Pop(k+1) is set to zero the proposed algorithm 
behaves as standard LMS algorithm. If the elements in tap weight vector is either zero or higher than 1/a, 
then there is no pull towards zero on any coefficients and the term Popi(k+1) is equal to zero. 
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2.2. Extraction of active and reactive load current component 

Figure 2 represents a detailed diagram of lo-LMS based adaptive filter. After each iteration, the 
active weight components are evaluated and updated. The load current, weight, and input signal are used to 
determine the error signal. The weights adjust themselves to filter out unwanted harmonics in the load current 
and provide fundamental component of reference weights. Flowchart for the extraction of components is as 
shown in Figure 3. The weight vector of active component of a load current of phase ‘a’ is calculated as, 


Wpalk + 1) = Wpalk) + Upael(k) x(k) + Pfo(Wpalk)) (10) 
e(k) = i(k) — Upa(k)Wpalk) (11) 


In the same manner, active weight vector of load currents in the phases ‘b’ and ‘c’ are calculated as, 


Wpp (k + 1) = Wyn (k) + Mppe(k) x(k) + p&o(Wpn(k)) (12) 
Wyck + 1) = Wyc(K) + Upce(k) x(k) + p&o (Wye (k)) (13) 
Vac" 
Wsp Vide 
7 Wpa (kK + 1) = Wya (kK) + Up (rm) x(n) + po (Wa (k)) 
Upa iba = Wsp * Upa 
ius Wpp (k + 1) = Wyp(k) + Upp er) x(2) + pFo(Wps (K)) 
Ue tgp = Wap * Ups | di 
Reference sz 
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Figure 2. Control structure of the proposed system 


Furthermore, the weight vector of reactive component of phase ‘a’ load current is estimated as, 
Waalk + 1) = Waalk) + Ugae(k) x(k) + pSo(Waalk)) (14) 
In the same manner, weight vector of reactive components of phases ‘b’ and ‘c’ load currents are calculated as, 
Wap (k + 1) = wan (k) + Mgel(k) x(k) + pFo(Wan(k)) (15) 
Wac(k + 1) = Wack) + Mgce (kK) x(k) + pSo(Wac(k)) (16) 


2.3. Extraction of active and reactive loss components 
Weight of active as well as reactive loss component is given as, 


Wep (kK + 1) = Wep(k) + Kp(vac(k + 1) — vac(k)) + kiVac(k + 1) (17) 


Weq (Kk + 1) = Woq(k) + kp (ve (k + 1) — ve(k)) + kive(k + 1) (18) 
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Figure 3. Flowchart of proposed algorithm 


2.4. Estimation of unit templates 
In phase as well as quadrature unit template can be calculated as, 


— Ya = 2b — Ve 
Upa = me Upp = a and Uj, = , (19) 
[2 
= 2 2 2 
Where 1%, = 5 (v2 + vp, + v2) (20) 
Upb—Upc V3Upa , UpptUpe —V3Upa , Upb—Upe 
Ugg = > Ugyp = —— + —" and ug, = —— + 21 
ae v3” a> Mere qb a toa Gy 


2.5. Generation of reference current 

The weight vectors are averaged to determine the equivalent weight for active and reactive component 
of current. The averaging of weight vector aids in minimizing any imbalance in the current components. The 
total active weight vector wr, and reactive weight wr, of grid current reference can be given as, 


Wrp = Wipa + Wep (22) 
Where Wp is total average active component of load current and w,, is active loss component, 


(23) 


Where W; qq is total average reactive component of load current and W,, is reactive loss component, 


3. RESULTS AND DISCUSSION 

The proposed system as well as adaptive LMS based algorithm is developed and simulated in 
MATLAB Simulink software using the simpower system (SPS) toolbox for simulation. The SPV array is 
designed for a maximum power output of 7.5 kW, and the PV inverter is rated at 7.5 KVA. The simulation of 
SPV array is as per [25] Also, the design of its various components is carried out in detail, such as value of 
the DC bus voltage and capacitor, the interface inductance, and the ripple filter. Under nonlinear operating 
conditions (balanced and unbalanced) and variable solar irradiance, the performance of the presented system 
is simulated. Various simulation parameters are as in Table 1. 


Table 1. Simulation parameters 


Parameter Value 
LMS covergence coefficient 1 5e3 
ZA strength and range a le? 
Smoothing parameter B 0.99 
Carrier frequency 10 KHz 
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3.1. Simulation results 
3.1.1. Performance of the system operating to support unity power factor (.) 

The proposed lo-LMS control algorithm is simulated for non-linear unbalanced load condition. The 
unbalanced condition is created by disconnecting load in phase C at 0.5s and its effects on source current, 
source voltage and compensator current waveforms are observed. As shown in Figures 4(a) and (b) the PV 
inverter is supplying extracted PV power to the connected load and or grid as well as acting as DSTATCOM 
compensating harmonics injected into the system by non-linear load and trying to maintain sinusoidal grid 
currents. The active and reactive both powers will be supplied from PV array. The DC link voltage and PCC 
voltages are also maintained constant. Figure 4(c) also shows various internal control signals. 
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(c) 
Figure 4. Performance of the system under unbalance non-linear load for maintaining power factor (a) Vabc, 


Iabc, load current (a and c phase), converter current (a and c phase), (b) Vac, Ipy and Py, and (c) intermediate 
control signals 


3.1.2. Performance of the system operating to support ZVR mode 
The inverter has the capacity to exchange both active and reactive power. It utilizes the remaining 
capacity of inverter after active power exchange. In addition, when PV power is not available (either due to 
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whether conditions or in nighttime), the PV inverter will now utilise as active filter to support required 
reactive power. Grid is supplying active power and grid current is maintain sinusoidal as PV inverter now 
manages reactive power. Figures 5(a) and (b) shows dynamic performance of the system when operating 
under ZVR mode. 


3.1.3. System dynamics with variable solar irradiance 

The behaviour of the system for variable solar insolation is as shown in Figure 5(c). At 0.15 s, the 
solar irradiance is increased from 500 W/m? to 800 W/m?. The maximum power is extracted from SPV array 
by MPPT algorithm and the extracted power is supplied to the load as well as grid. After increase in PV power 
the grid current increases, slightly as more power is now injected into the grid. Load is maintained constant so 
surplus power is injected into the grid and active grid power is negative indicating reverse power flow into the 
grid. The harmonic spectrum of grid current and load current is compared in Figures 6(a) and (b). 


3 509 ee Se ee ee ae ee ee ee a ee 
>° - 


times) 


Figure 5. System performance under various operating modes (a) dynamic performance under ZVR mode, 
(b) Active power distribution under ZVR mode, and (c) system performance under variable PV power 
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Figure 6. Harmonic spectrum (a) source current and (b) load current with non-linear load 


3.1.4. Performance comparison of proposed algorithm with dq and LMS algorithm 

The proposed ly-LMS control algorithm is compared with dq and LMS algorithm as shown in Figure 7. It 
is observed that the settling time is less for lo-LMS algorithm as compared to other algorithms under consideration. 
The value of the fundamental active component extracted during unbalanced load conditions by lo-LMS is more 
accurate than dq and LMS algorithm for same condition. Total harmonic distortion (THD) for all three algorithms 
is as shown in Table 2. THD is 3.8%, 4.1% and 3.37% for dq, LMS and lo-LMS algorithms respectively. 
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Figure 7. Performance comparison of dq, LMS and lb-LMS 


Table 2. Comparision of various parameters for dq, LMS and lo-LMS algorithm 


Performance parameter DQ LMS lo-LMS 
THD 3.8% 4.1% 3.37% 
Complexity Less Medium Medium 
Settling time 0.18s 0.17s 0.1s 
Accuracy Good Good Better 
No. of computation Medium Less Less 
Oscillation Medium High Slightly less 


3.2. Laboratory setup 

A prototype of the proposed system is developed in the laboratory, and it is tested for various 
conditions is as shown in Figure 8. An autotransformer is used to step down 400 V, 50 Hz, three-phase grid 
voltage to 110 V. A 6 mH inductor is connected to 110 V for realising the transmission line. Six PV arrays 
(WAREE make, each of 150 Wp) are connected in series for 900 Wp output. Switching frequency ripples are 
reduced by an RC filter connected after the three-phase H bridge inverter. Hall effect voltage/current sensors 
(NITECH made) are used for sensing various quantities like DC voltage, grid voltage, and grid/load/PV 
currents [23]. The proposed adaptive lo-LMS for reference current generation is implemented using ARM- 
CORTEX M4 microcontroller (STM32F407VG). For recording the performances of the proposed system, a 
four-channel DSO (Agilent make) is used. The system is tested under various conditions and the results are 
presented as below. The hardware results for various conditions are as shown in Figure 9. The system is first 
tested for only active power transfer from the PV array to the grid. Figure 9(a) shows inveter currents for 
different solar irradiance and observes that inverter currents change accordingly. Figure 9(b) shows a 
dynamic change in grid current for a change in solar irradiance as PV power is directly supplying load/grid 
and also DC link voltage is maintained constant with PV power variation (PV power variation is achieved by 
covering the PV array with a semi transparent sheet). Only active filtering mode is shown in Figure 9(c) and 
Figure 9(d). Inverter current follows the reference current generated using proposed algorithm in Figure 9(c) 
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and in 9(d); its effect on grid current is shown. As is clear that after inverter current follows reference, the 
grid current becomes as PV inverter supplies required reactive power. Figures 9(e)-(f) shows hybrid mode in 
which active and reactive both powers are supplied by the PV inverter. PV power variation for the conditions 
like low load-high PV power as well as high load-low PV power is depicted in Figure 9(e). Figure 9(f) shows 
dynamic condition of all modes i.e., upto region 1, only real power is injected into grid, 2" region shows 


condition where both real and reactive power is supplied by PV inverter and at 3" region PV aray is removed 
so now inverter will act as active filter only. 


2. Resistive load 
3. Variac 

4. Rectifier 

5, Line inductor 


6. VSC 
1 7. Sensor card 
may a2) °P 8. STM32F407VG 


9. DSO 


510. SPV array 
~ 


| 1, SPV array termination 
h 


Figure 8. Hardware setup for proposed algorithm 
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Figure 9. Experimental results under different modes (a) steady state performance for PV mode, (b) dynamic 
performance for PV mode, (c) active filter mode (inverter current), (d) active filter mode (grid current), 
(e) hybrid mode (steady state performance), and (f) hybrid mode (dynamic performance) 
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4. CONCLUSION 

The proposed adaptive lo-LMS algorithm is implemented in the PV-DSTATCOM system for 
mitigation of various PQ issues. The proposed algorithm has been tested for conditions of load ubalancing 
and variable solar irradiance and the results are found satisfactory. The adaptive lo-LMS algorithm has the 
ability to adapt to time-varying signals as it generates proper control signals under varying load and 
atmospheric conditions. Zero attractor is added to the update equation of tap weight vector, which increases 
the convergence speed. From the results, it is also observed that the PQ of the system is maintained while it 
continues to supply maximum power extracted from PV array. Results are verified on low cost laboratory 
prototype developed using STM32F407VG microcontroller. 
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